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ABSTRACT 


A numerical simulation model has been developed for the prediction of 
fallout from subsurface nuclear detonations that produce craters through 
spall and the action of the cavity gas. The physical processes modeled are 
atmospheric transport, lateral eddy diffusion, and gravitational sedimenta- 
tion of radioactive particulates. This "cratering fallout model" is normalized 
to the observed external gamma dose rate fields of the Sedan (100-kt) and 
Danny Boy (0.43-kt) cratering events conducted at the Nevada Test Site. Cal- 
culations of the fallout patterns for additional events, used for testing the 
prediction capability of the cratering fallout model, indicate that the model 
gives estimates of the external gamma dose rate at H+ 1 hours with a maxi- 
mum error of a factor of 2 to 3 in the gamma dose rate versus distance along 


the hot line of the pattern. 


I. INTRODUCTION 


During the past few years, a small but continuing effort has been ex- 
pended in developing a model for predicting fallout from subsurface nuclear 
detonations. In this report, we (a) summarize the experience in predicting 
fallout from surface bursts that is transferable to the construction of a fallout 


model for subsurface detonation, (b) describe the development of the fallout 


model for subsurface nuclear detonations (the cratering fallout model), 

(c) illustrate by means of an independent test case the predictive capability 
of the model, and (d) discuss some of the problems associated with the pre- 
diction of fallout from row-charge subsurface nuclear detonations. This 
paper is limited to research or development in which the investigator has 
been personally involved within the Plowshare Program of the Lawrence 


Radiation Laboratory. 


Il. THE CRATERING FALLOUT MODEL 
A. Basis of Model Construction 


In order to predict the area affected by radioactive fallout from a supe 
surface detonation of a nuclear explosive, and to predict the gross external 
gamma dose rate in the surface fallout pattern, we need knowledge of the 
following: 

1. The height of the base and top, and the radius, of each radioactive 
cloud (i.e. , the main cloud and base surge cloud) formed by the detonation at 
the time the clouds cease to rise inthe atmosphere. (This time is defined 
as the time of cloud stabilization. ) The cloud heights are prescribed in terms 
of height above surface zero. 


2. The total yield of the explosive, W the fission yield of the ex- 


7? 


plosive, W,, the depth of burial of the explosive, z, the fraction of the fission- 


f’ 
product gamma emitters expected to appear in the fallout pattern beyond the 
(estimated) radius of direct ejecta, Bo and the equivalent fission yield 
gamma needed to simulate gamma dose from induced activities. 


3. The activity-particle size distribution in the main cloud and the 


base surge, and the fraction of By in each of these clouds. 
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4, The terminal fall velocity of the fallout particles (in still air) as a 
function of particle size and height in the atmosphere. 

5. Ina fallout calculation involving two clouds, main cloud and base 
surge, the calculation of the fallout pattern for each cloud is done separately. 
The total surface fallout pattern of the event is found by summing the patterns 
from the base surge and main clouds. Thus, the time and space prediction of 
the horizontal wind at the level of each cloud top is required, along with the 
specification of the wind-shear tangential and normal to the wind for the layer 
through which particles fall. 

6. The effect of horizontal eddy diffusion on the growth of the horizontal 
radius of the disks of radioactive particles as the disks fall earthward. (The 
initial debris cloud is subdivided into disks of debris in the model as a func- 
tion of initial height in the cloud and particle size. ) 

The above information provides the simplest but still adequate basis 
for constructing a cratering fallout model. In addition to this information, 
we assume that the fission. product radioactivity is unfractionated and that 
1 kiloton of unfractionated fission products spread uniformly over 1 square 
mile corresponds to an (H + 1)-hour dose rate of 3380 R/hr at a height of 3 
feet above an infinite plane [ Higgins, 1963]. This normalization constant of 
3380 R/hr is corrected for terrain shielding by a factor of 0.8. 

Certain atmospheric processes or effects have been neglected in de- 
veloping the cratering fallout model, just as these same effects were neglected 
in developing the land surface-burst fallout model [Knox, 1964a]}: 

7. The effect of synoptic-scale vertical motions on the vertical dis- 
placement of the falling disk of particles. 

8. The effect of the mean divergence of the horizontal wind on the radius 


of the disk of falling particles during descent. 
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9, The effect of vertical eddy diffusion. 
10. The time from detonation to cloud stabilization. 
_ 11. The development of subsynoptic-scale wind systems. 
12. Changes in time and space of the shear tangential and normal to 
the horizontal wind at cloud top level. 
13. The effect of water products of condensation on the size and fall 
rate of radioactive particles. 


We consider now, more fully, items one through six. 


B. Fraction of Gamma Activity in the Close-In Pattern 


In order to calculate the close-in fallout pattern from a subsurface detona- 
tion, the fraction of the gamma activity produced by the fission yield and 
appearing in the close-in fallout pattern F. must be known. To date, ex- 
perimental data from Sedan, Teapot ESS, Jangle U, Neptune, and Jangle S 
(and other surface burst events) have provided a basis for estimating F\.. 
Integration of the measured fallout patterns for these events has been per- 
formed [Williamson, private communication, 1964] from the radius of direct 


throwout to the limit of the measured pattern, with results as shown in Table I. 


Table Il. Experimental FB data. 


Event W, kt z, ft By Medium 
Sedan 100 635 ~0.10 _ Alluvium (A) 
Teapot ESS 1,2 67 0.46 A 
Jangle U 1.2 17 0.64 A 
Neptune 0.115 100 0.005 Tuff 
Jangle S 1.2 0 0.50 A 
Danny Boy 0.43 109 0.04 Basalt 


Blanca 19 835 0.0005 Tuff 


ev 


These FB. calculations were performed assuming a normalization constant 
of 3380 R/hr per kt of fission yield per square mile at (H + 1) hours and a 
terrain shielding factor of 0.8. The graphical presentation of F. as a func- 
tion of 2pw!3 is shown in Fig. 1. In constructing the experimental BR 
curve, we have been guided by: 

1. An asymptote of FB. = 0.50 for z=0, suggested by Jangle S data in 
the Table I and supported by previous work [Nordyke and Wray, 1964]. 


1/3 ~ 330 tt/ut!/3 


2. An asymptote of no dynamic venting for z/W 
supported by experimental evidence [Higgins, private communication, 1964] 
and by studies of containment physics [Knox, unpublished Internal LRL 
Memorandum, 1963]. 

It should be noted that the value of BP. (estimated by pattern integra- 
tion) for the Sedan event departs considerably from the curve fitted to the 
whole sample of BF. data. This departure could arise because of an B 
yield dependency that is inadequately known at this time. In the absence of 
knowledge of such an BF dependency, re estimates made from the By curve 
in Fig. 1 for high-yield cratering events (of the order of 100 kt and above) 
should be considered as uncertain by a factor of 2. It should be further 
stated that F. estimates derived from Fig. 1 tacitly assume that the nuclear 
explosive is fully tamped and that no additives for fission-product gamma 
radiation suppression have been placed around the explosive. 

It should also be mentioned that the maximum value of Be (0.75 at 
2/ wil? = 30) occurs at very nearly the same value of 2/w/ 3 as the maxi- 
mum base surge radius (crosswind) in a neutral atmosphere, reported earlier 


[Knox and Rohrer, 1963]. The author believes that this coincidence of 


maxima at z/w'!3 = 30 is physically consistent. 
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Fig, 1. Fraction of the gamma activity appearing in the close-in fallout 
pattern, 
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C. Cloud Geometries 


The geometrical definition of the cloud top, base, and radius for the 
main cloud and base surge cloud at the time of cloud stabilization is shown 
in Fig. 2. In addition to the definitions given in Fig. 2, the height of the base 
of the main cloud is defined as being equal to the height of the base surge 
cloud top Hi, in the model. 

Evidence suggests that the geometry of these two clouds at the time of 
stabilization is a function of the total explosive yield, the material in which 
the detonation occurs, depth of burial of the explosive, and the meteorological 
conditions existing during the development of the clouds. [Knox and Rohrer, 
1963]. At present, the cloud geometry parameters (Ry, HH Ron H,,) must 
be evaluated experimentally as functions of total yield and depth of burial. 
Reasonable samples of experimental data exist for alluvium and basalt mate- 
rials, An example of one of the most useful summaries of cloud geometry 


data for alluvium is shown in Fig. 3b-g- [W. Day, U. 5. Army Corps of 


Engineers Nuclear Cratering Group, LRL, private communication, 1964]. 


Day's summary utilizes all of the known cloud geometry data from high ex- 


plosive and nuclear explosive detonations conducted by the U. S. Atomic 


-Energy Commission in alluvium and basalt. Implicit in Day's summary is the 


assumption that a high explosive and nuclear explosive detonated in the same 
material at identical depths of burial and under similar meteorological con- 


ditions produce the same cloud geometries. 


* Figure 3a is a computational aid to the acquisition of input to Fig. 3b-g. 
In these figures, DOB denotes depth of burial, and D, denotes the depth of 


apparent crater. 
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Fig. 2. Definition of cloud dimensions and symbols. 
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Fig. 3a. Computational aid to the acquisition of input data to the cloud dimension graphs, Fig. 3b-g. 
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Fig. 3b. Base surge radius as a function of yield (total) and the parameter DOB/ dD, (alluvium). 
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Fig. 3d. Base surge height as a function of the yield (total) and the parameter DOB/ D, (alluvium). 
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Fig. 3f. Main cloud radius as a function of the yield (total) and the parameter DOB/ dD, (alluvium). 
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Fig. 3g. Main cloud height as a function of the yield (total) and the parameter DOB/ dD. (alluvium). 
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D.. Activity-Particle Size Distributions 


ry) 


In a typical subsurface nuclear detonation in alluvium, two clouds are 


formed. The main cloud is composed mostly of vented cavity gas and particu- 


lates (originating either from condensation or injection and entrainment of 
soil).. The base surge cloud is composed of ejecta and suspended fine par- 
ticulates. For a nuclear cratering shot in alluvium, it is assumed that 80% 
of the F gamma activity is in the main cloud and 20% in the base surge 
cloud. The 0.8 Be main cloud activity is assumed to be subdivided between 
two log-normal activity-particle size distributions. The first activity- 
particle size distribution contains the activity 0.8w (LF, and is characterized 
by the mean ln f(1) and the standard deviation O61). The second activity - 
particle size distribution contains the activity 0.2 Wy l2F and is charac- 
terized by the mean In rl?) and the standard deviation o, (2). The activity 
of the first distribution is assumed to be hemeaensously mixed through the 
whole main cloud, while that of the second distribution is assumed to be 
homogeneously mixed in the lower fifth of the cloud. A similar prescription 
of activity vs particle size is used for the base surge cloud. Figure 4 sum- 
marizes the parameters governing the activity-particle size distributions in 
the two clouds at the time of cloud stabilization. These parameters govern- 
ing the activity-particle size distributions in the cratering fallout model have 
been determined by mathematical experimentation with the model in the re- 
creation of the observed fallout patterns for the Sedan and Danny Boy events. 


Results of these calibration calculations will be discussed in a later section. 


E, Terminal Fall Velocity of Fallout Particles 


The vertical fall velocities of the fallout particles are modeled as the we 


terminal fall velocities of smooth spheres of density 2.5 g/cc inan LC. A.O, 


0.8 Wm(l)Fe 


Wm (1) 
Ln Fp (1) 


om (1) 


Wm (2) 
wp (2) Lnim(2) 
Rn rp (2) em(2) 
a (2) : = 
0.2 wp (2) Fo Se a 


SURFACE 
ZERO 


Fig. 4. Schematic drawing of an idealized cloud from a subsurface detonation showing spatial relation- 
ships of the activity-particle size distribution assumed in the model, 
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standard atmosphere as computed by McDonald [1960] for both the Stokes 

law region and the aerodynamic region (wherein the Reynolds number exceeds 
1). If in reality some fallout particles are a cluster of small spheres attached 
to a large central particle, these complex particles are assumed to fall with 
the speed of the equivalent smooth spherical particle of the same mass. Ex- 
perimental evidence has been obtained by Rapp and Sartor [1957] to support 


this assumption. 


F, Specification of the Horizontal Wind Field for the Model 


The horizontal wind field that transports the debris disk centroids dur- 
ing their fall to the earth's surface may be specified in two ways in the crater- 


ing fallout model. 


Option 1. Idealized Wind Hodograph 


If H denotes height above surface zero, then the horizontal wind ea 


at height H for a simple wind hodograph (see Fig. 5) is 


: Alp) ; 
vy) = v yl.) go> - SH, - Hg 


where 
vf) can be specified as a steady state wind at the cloud top 
level A or v fA) can be specified as a function of 
time and space in either wind component form or by 
means of a stream function; 

A(p) characterizes the wind shear component tangential to the 
horizontal wind at cloud top level and is evaluated from 
shot time winds (in a diagnostic fallout calculation) or 
from preshot wind information (for a predictive fallout 
calculation), and it is held constant in time for the period 


~ of fallout deposition; 


ag 


~ 


Fig. 5. 
(Option 1). 
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Alp,) is the value of A(p) at cloud top level (normally it is set 
equal to 1); 

S is the wind shear conmouent normal to the horizontal 
wind at cloud top level, evaluated from shot time winds 
or preshot wind information (depending on the purpose of 
the fallout calculation) and held constant in time; 

nis the unit vector normal to v fA) in a right-handed 
system; 


p is the atmospheric pressure corresponding to H. 


Option 2, Arbitrary Hodograph 


In Option 2 the horizontal wind can be specified in wind component form 


for as fine a vertical interval as desired or for which wind information exists. 


G. Debris Disk Radius as a Function of Time 


To estimate the radius of a debris disk expanding by horizontal eddy 
diffusion during its fall to earth, we propose that the disk radius as a function 


of time R(t) be represented by 


‘ 11/2 
R(t) = (R, 4 + 2Dt') 


where t! is the distance traveled by the disk centroid divided by the mean 
horizontal wind speed in the layer through which the disk has settled, Ro is 
the debris disk radius at time of stabilization, and D is the horizontal eddy 
diffusion coefficient [Knox, 1964a]. The diffusion coefficient D is estimated 
as the Richardson's diffusion coefficient 0.2 x (2)4/8 where 7 is the standard 
deviation of the position of the particles (in the disk) from the disk centroid. 
Since Z is poorly known in nuclear debris clouds, it is usually set equal to 


R or 2k, if accelerated relative diffusion is to be approximated in the 


e,0 0 
fallout model. 


a 
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H, Physical Processes Simulated in the Model 


The principal physical or meteorological processes simulated by the 
cratering fallout model are: (a) the transport of the debris disks by the mean 
wind in the layer through which the disks are falling, (b) the relative advection 
of the debris disks by the horizontal wind field containing both speed and 
directional shear, and (c) lateral eddy diffusion which expands the disks falling 
earthward. 

The first two processes are simulated by calculating the transport of 
the disk centroids by the ambient horizontal wind field during the disk's 


earthward fall until the disk centroid is on the ground surface, using either 


‘Option 1 or 2 for the horizontal wind specification. The result of this calcula- 


(ion is the position vector from surface zero to the predicted ground position 
of each disk centroid and the time of arrival at the ground surface of each 
disk. Process (c), above, is modeled by the expression for R(t) which 
estimates the debris disk radius at the time of arrival of the disk centroid on 
the ground surface. | 

The debris disks tracked earthward are defined as follows: In each 
cloud (main cloud and base surge cloud) 11 disks of particle size ry are de- 


fined at each isobaric surface Po such that 


Inr,=Inr 4 Se 


I L410 ier! 


where 
b= 0, 1x 246 2% 5 105, 


rr is the radius of the largest particle modeled in the cloud, 

Tg is the radius of the smallest particle modeled in the cloud, 
Z. P = 

p PE aq (Pee Pe and: FP =4O9 Te. Bi Seay Oy 

Pp is the pressure at the base of the cloud at the time of cloud 

stabilization, 
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Pn is the pressure at the top of the cloud at the time of cloud 


stabilization. 


The (H +1)-hour external gamma dose rate (for a height 3 ft above an 


infinite plane) is calculated using the method of Batten, Inglehart, and Rapp 


[1960] modified to account for the effect of normal shear S and lateral eddy 


diffusive disk growth. 


The following quantities are calculated in the model and are output in 


the indicated modes 2 


(a) 
(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


Printout CRT display 
Position of surface zero x 


Dose rate (H + 1) at the predicted x 
ground position of each disk 
centroid for each cloud 


The predicted ground position . x x 
of each disk centroid for each 


cloud 


The envelope containing the area x 
affected by the fallout from each 


cloud 


The time of beginning and end of x 
fallout deposition at each ground 


position centroid for each cloud 


The isopleths of the (H + 1)-hour x 
dose rate for any specified interval 
of dose rate (each contributing cloud 


and total pattern) 


The (H + 1)-hour dose rate as a x 
function of distance along the hot 


line (each cloud and total pattern) 


*For a two-cloud nuclear cratering fallout problem about 0.25 min of 7094 


time is required for the calculations, and 2.4 min of LARC time is required 


for the CRT displays indicated above. 


«ed 
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Ill. DIAGNOSTIC CALCULATION FOR SEDAN 


In order to calibrate the cratering fallout model on the Sedan event, 
the observed shot time winds, the observed cloud geometry of the main cloud 
and base surge cloud, the estimated BR = 10%, and the appropriate fission 
yield were input to the model. A first '' guess" of the 12 activity-particle 
size parameters, discussed in Section IV, was also input to the model. The 
main cloud was assumed to contain 80% of Fo the base surge cloud 20%. By 
mathematical experimentation, the parameters governing the activity-particle 
size distributions have been determined from the observed Sedan gross gamma 
fallout pattern. The values of these parameters are as follows: 
w (1) = wiQ1) = 0.9 
w, (2) = wy (2) a Oe 
be Ch) dnee (e209 
In r (2) = In r,(2) = 5.0 
o(1).= o,(1) = 0.69 
o,(2) = o,(2) = 0.59 
Figure 6 shows the calculated and the observed gamma dose rate (at H + 1) 


from fission products vs distance along the hot line for the Sedan event. 


IV. DIAGNOSTIC CALCULATION FOR DANNY BOY 


The Danny Boy event was a 0.42-kt nuclear cratering detonation em- 
placed at a depth of 109 ft in dry basalt. The observed gamma fallout pattern 
for Danny Boy has been published [Nordyke and Wray, 1964]. By mathemati- 


cal experimentation with the cratering fallout model, we can adjust the 


“The observed (H + 1)-hour gamma dose rate for fission products derived 
from the observed total gamma dose rate at H+ 24 hours, assuming that 52% 


of the (H + 24)-hour gamma dose rate was due to tungsten [Knox, 1964b]. 
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activity-particle size distribution parameters in order to duplicate the ob- 
served Danny Boy fallout pattern. In this event no visible main cloud was 
observed [Nordyke and Wray, 1964]. Thus in the diagnostic calculation with 
the cratering fallout model we assume that 100% of the FB. gamma activity is 
in the base surge cloud, and a preliminary value of F. = 0.05 is assumed. 
The activity-particle size distribution parameters determined in this calcula- 


tion are given below: 


wp (1) = 0.9 
wy, (2) = 0.1 
In rp) = 3.0 
In vr, (2) = 5.7 
o,,(1) = 0.69 
o,,(2) = 0.59 


Figure 7 shows the calculated (H + 1) dose-rate-vs-distance curve from the 
model, and the observed (H + 1) dose rate vs distance. As was previously 
mentioned, the value of RB used in the cratering fallout model calculation 
was 0.05; if the value of F, = 0.04 (reported later, Nordyke and Wray [1964]) 
had been used, the agreement between the calculated (H + 1) dose rate and 
the observed (H + 1) dose rate vs: distance would have been better than shown 
in Fig. 7. 

Recently, effort has been expended to obtain machine capability of 
plotting the calculated (H + 1) dose rate patterns for an arbitrarily selected 
interval of dose rate. Figure 8 shows the machine-plotted dose rate pattern 
for the. Danny Boy diagnostic calculation. The computer pattern breadth at 
25 miles downwind is about 7 miles, and the breadth of the observed pattern 
at this distance is 5.5 miles. It should be noted that the closure of the isodose 


rate line of 1074 R/hr at (H + 1) in the plotted pattern is artificial and is the 
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Fig. 8. The machine-plotted (H + 1)-hour gamma dose rate pattern for the 
Danny Boy event (diagnostic calculation). 


result of the logic used for the computer plotting, rather than the result of 


the cratering fallout model. All isodose rate lines will apparently be closed 
at the downwind edge of the pattern if computed dose rate information is in- 
sufficient for their appropriate extension downwind. The observed (H+ 1) 
gamma dose rate pattern for Danny Boy is shown in Fig. 9 for comparison 


with the computer -plotted pattern, 


V. INDEPENDENT TEST OF THE MODEL FOR TEAPOT ESS 


An independent test calculation of the cratering fallout model (calibrated 
on Sedan) was performed using the Teapot ESS event. The observed winds 
at shot time, the observed cloud geometry, the published fission yield of 1.2 
kt, and F. = 0.85 (an early preliminary value of F.) were input to the model. 
Figure 10 shows the calculated and observed (H + 1) gamma dose rate as a 
function of distance along the hot line of the pattern. A comparison of these 
two dose-rate-vs-distance curves shows that the largest error between cal- 
culation and observation is of the order of a factor of 2.5. An examination 
of the radiosonde observation (NTS) near shot time indicates that the vertical 
temperature distribution of the layer through which the particles were falling 
was Slightly superadiabatic. Under such conditions, it is possible that the 
vertical eddy diffusion on the Teapot ESS shot day was larger than average. 

If enhanced vertical eddy diffusion were operative, the dose rate would be 
slightly less near ground zero and enhanced downwind. 

Other model confirmation tests have been performed. These include 
calculations with the surface burst version of the model on the Apple II and 
Zucchini events [Knox, 1964q and several other atmospheric events whose 
results Pewiain classified at this date. The results of these additional tests we 


were satisfactory. 


af 
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PRELIMINARY ISODOSE CONTOURS 
DANNY BOY EVENT 


Readings in R/hr-H+lhr 
Data taken from: 


NDL Ground surveys (close in) 
EGG Aerial survey (intermediate range) 
USGS Aerial survey (long range) 


GLL -625-16926A 


Fig. 9, Preliminary gamma isodose rate contours observed at H+ 1 hours, 


from Nordyke and Wray [1964]. 
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Fig. 10. The calculated and observed gamma dose rate at H+ 1 hours asa 
function of distance along the hot line of the pattern (predictive calculation), 


Teapot ESS. we 
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VI. COMMENTS ON THE PREDICTION OF FALLOUT FROM 


ROW CHARGE EVENTS 


In principle the fallout pattern from a row of subsurface nuclear detona- 
tions may be estimated with the cratering fallout model, provided the model 
input parameters are adequately known for row charge events. In the current 
state of knowledge, the cloud geometry input parameters and the value of B, 
appear to be the most difficult to specify. Study of time-lapse photography 
of the main cloud and base surge cloud evolution from past H. E. row-charge 
tests give significant experimental information on the cloud geometry param- 
eters. For example, empirical methods of predicting crosswind radius and 
height of the base surge originating from H. E. row-charge events have been 
discussed [Knox and Rohrer, 1963]. In this latter study, it was shown that for 
five equal-weight charges, equally spaced and emplaced at the same depth of 
burial, the resulting base surge cloud has approximately the same radius as 
the base surge cloud from a large single shot of the same total yield and at 
the same scaled depth as the smaller charges. The base surge height for a 
59-charge row event is approximated, reasonably well, by scaling the height 
of surge cloud for a single-charge event by the 0.2 power of the total yield of 
the row charge events. 

A first approximation to the geometries of the individual main cloud may 
be obtained by treating each main cloud independently and estimating the cloud 
top and radius of the individual cloud from the work of Day [W. Day, U. S. 
Army Corps of Engineers Nuclear Cratering Group, LRL, private communica- 
tion, 1964]. Examination of the H. E. row-charge documentary photographs 
of the Rowboat, Dugout, and pre-Buggy events indicates that such an approxi- 
mation is reasonable. This approximation can, of course, be in error if 


main cloud interactions occur. There is a need to evaluate the uncertainty in 


waoe 


fallout prediction for multiple-charge events in cases where cloud interactions 
‘lead to the injection of radioactivity at levels higher in the atmosphere than 
predicted. | 

Concerning the BR, for row charge events, preliminary experimental 
results for H. E. single and row-charge events have been reported previously 
[Graves, Wray, and Pierce, 1963]. In this study, results of the measured 
vented fraction of lia’ *” tracer from single and row-charge H. E. events were 
given. The experimental evidence biapeste that the vented fraction from row- 
charge events may be about twice that for single-charge events. One of the 
more pressing needs for fallout prediction from row-charge events is the 
establishment, either through experiment or theory, of the dependency of the 
Fy on yield, depth of burial, or charge spacing. 

However, the development of computer aids for the row-charge fallout 
prediction problem can proceed independently of the solution of the two above- 
cited problems. To this end, the capability of plotting fallout patterns for 
multicloud and/or multidetonation events has been developed Figures 11 and 
12 show, respectively, the (H + 1) dose rate patterns computed for 10 Danny 
Boy detonations on an east-west line with charge centers separated by 33.5 
meters, and 10 Danny Boy detonations on a north-south line with charge 
centers separated by 33.5 meters. For the purpose of these calculations, 
each detonation is assumed to vent 5% (e. g., By = 0.05), the individual clouds 
are assumed to be that of Danny Boy, and the input wind for each cloud fallout 
problem is assumed to be the shot time wind for Danny Boy. These two row- 
charge fallout patterns, although detectably different as determined from the 
printout, appear to be very similar. These calculational results suggest 
that the fallout pattern from a small-scale, row-charge event is not sensitive 


to the orientation of the wind to the alignment of the charges. 


we 


204 


‘ ’ 
a a Ce 


10 DEVICES IN A HORIZONTAL LINE 33.5 METERS APART 
Dae) LNPRisz LNPGOSE-7 Sis-69 =-59 FORTRAN CODE 
eee CMM OLE TC 6 CRIES 00 0 MMEDIUILG,® 0 0 GLL~6410- 3068 


“ya 72958 OUIAIN «owi-THITTE C2 YMAX = 6 -S0GG5E OH 


Fig. 11. The (H + 1)-hour gamma dose rate pattern computed for 10 Danny 
Boy detonations on an east-west line with charge centers separated bv 33.5 
meters, 


-34- 


2 Z 
eee eneesenesetan 


: - : ¢ . . 
eee ee eo ee ee ) Ce) eae nah aeae eevee eeeeeonetaosenstoesreseesess 
. . a . 


age ie he ee ee a oe ee TO ek Se NU eB ES a a Co cree ee ie ea ee Me ene ote el gO ag tm anh 
Pr er eo ee re CTC ee ee ce i a a ee a er 


ee ae ee ee ee 


baer act TODEVIGES IN A VERTICAL LINE 33 SPETERS APART Bee ee é 
DOSE RATE CO} OOO OO OT ~~» GLL-64.10-3069 |." 
TEN BILES BETK- GRID LINES ee 
TAIN A OMEOEE OF RAL S“1TNS7E GUTHIN . N-RONSTE 03 YRAK . 6-9368NE OF 


Fig. 12. The (H+ 1)-hour gamma dose.rate pattern computed for 10 Danny 
Boy detonations on a north-south line with charge centers separated by 33.5 
‘meters, 
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VII. CONCLUSIONS 


In conelusion, it may be stated that the cratering fallout model developed 
gives reasonable results for the prediction of the area of the pattern, geometry 
of the dose rate contours, and (H+ 1) dose rate vs distance along the hot line. 
The independent testing of the model should be extended to include more than 
the Teapot ESS case shown and the several others cited. 

Meteorological improvements in the cratering fallout model could well 
include the following items: (a) the prediction of synoptic-scale changes in 
the normal and tangential shear components, (b) better understanding of the 
growth of the debris disks through horizontal eddy diffusion during their earth- 
ward fall, (c) the prediction of diurnal change of wind at low levels and close 
to ground zero, (d) the inclusion of the effects of topography on the evolution 
of the horizontal wind field, (e) the effect of terrain-induced circulations on 
fallout deposition, and (f) a better solution to the cloud rise problem for 
cratering detonations. 

However, before embarking on the development of the above-cited 
meteorological improvements in fallout prediction, one must consider potential 
improvements which may come from other areas. Areas of promise are 
cratering physics and filtration theory for the vented-fraction problem, special 
emplacement for control of vented fraction, and the improvement of nuclear 


explosives, 
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